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ABSTRACT 


In  this  thesis,  rational  interpolation  operators  are  con¬ 
structed  from  polynomial  interpolation  operators  in  such  a  way  that  the 
rate  of  convergence  is  improved.  A  short  description  of  the  method  to 
be  used  is  contained  in  Chapter  1. 

Chapters  2  and  3  deal  with  the  background  material  required 
to  apply  the  method.  A  Korovkin-type  theorem  which  appears  in  Chapter 
2  forms  the  means  to  determine  rates  of  convergence  for  the  rational 
interpolation  processes.  The  required  information  regarding  the  poly¬ 
nomial  interpolation  processes  to  be  used  is  in  Chapter  3. 

The  new  interpolation  processes  are  constructed  in  Chapters 
4  and  5,  and  their  rates  of  convergence  determined.  The  results  in 
Chapter  4  are  of  a  more  general  nature,  while  those  of  Chapter  5  yield 
a  better  degree  of  approximation. 

Chapter  6  contains  some  suggestions  for  ways  that  these 
results  could  be  extended  or  improved. 
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CHAPTER  1 


INTRODUCTION 


The  problem  of  finding  interpolation  processes  which  provide 
good  tools  for  approximation  is  one  which  has  been  well  studied,  and  a 
number  of  such  processes  have  been  found.  It  is  natural  to  attempt  to 
modify  known  processes  in  order  to  improve  on  the  degree  of  approxi¬ 
mation  obtained.  This  thesis  contains  a  number  of  such  attempts.  The 
general  method  of  modification  is  natural  enough,  and  could  likely  be 
applied  to  interpolation  processes  other  than  those  considered  in  this 
thesis : 

Let  be  a  positive  linear  interpolation  operator  defined 

on  C [ —1 , 1 ]  by 

V(X)  =  X  8nk(x)f(Xnk)  > 

k=l 

where  -l<x  <  • • •  <  x  ,  <  1  and  the  g  .  are  nonnegative  con- 
—  nn  nl  ~  nk 

tinuous  functions  defined  on  [-1,1].  Suppose  that 

I  snk(x>  =  1  > 

k=l  n 

so  that  H  1  =  1.  Then  the  degree  to  which  H  approximates  continuous 
n  n 

functions  is  largely  dependent  upon  how  rapidly  |g  ^(x)|  decreases 
and  |x-x  increases.  This  dependence  is  expressed  precisely  by 
Corollary  2.1. 

Presented  with  such  an  operator  H^,  where  | gn^ Cx) |  0 

as  |x-x  I  grows,  one  can  often  realize  an  improvement  by  considering 
nk 

-  1  - 


' 


2. 


instead  an  operator  L  : 

n 


n 


L  f(x)  = 
n 


I  mnk(x)f(xnk) 
k=l 

£  mnk(x) 


k=l 


2  2 

where  m^Cx)  =  g^^(x),  or  something  very  like  g^^(x).  Then  m^(x) 

will  approach  zero  more  rapidly  than  g  .  (x)  as  x  moves  away  from 

nk 

x^.  Thus,  will  be  superior  to  as  a  tool  of  approximation, 

provided  that  the  minimum  value  achieved  by  the  denominator  does  not 
approach  zero  too  quickly  (or,  in  the  best  case,  is  bounded  away  from 


zero).  does  represent  an  increase  in  complexity  over  -  e.g., 

if  is  a  polynomial  operator  then  the  images  of  are  rational 

functions.  However,  the  denominator  of  L^f  is  independent  of  f,  and 
so  need  be  calculated  only  once  for  each  n. 


This  method  will  be  applied  to  operators  connected  with 
Hermite-Fej er  and  related  interpolation  processes,  which  will  be 
discussed  in  Chapter  3.  A  basic  paper  in  this  field  of  study  is 
Fejer  [5].  The  precise  effectiveness  of  the  H-F  process  was  discovered 
by  R.  Bojanic,  and  is  given  by  DeVore  [4,  p.232].  Questions  concerning 
slight  modifications  of  the  process  were  considered  by  Turdn  [11]. 

Other  modifications  and  extensions  were  defined  and  studied  in  Balazs 
and  Turan  [1,2,3];  Meir,  Sharma,  and  Tzimbalario  [6];  Prasad  and  Varma 
[7],  and  Sharma  and  Tzimbalario  [8]. 

The  nodes  of  the  processes  considered  will  be  the  zeros  of 
the  classical  ultraspherical  polynomials.  Information  about  these  can 
be  found  in  Szego  [10],  and  will  be  given  where  it  is  required. 
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CHAPTER  2 


A  KOROVKIN-TYPE  THEOREM 


In  this  chapter  will  be  stated  and  proven  a  Korovkin-type 
theorem  due  to  Mond  and  Shisha,  and  then  the  result  will  be  restated  in 
a  form  in  which  it  can  be  applied  to  the  operators  to  be  considered. 

Definition  2.1:  Let  L  be  a  linear  operator  defined  on  a  set 
D  of  real-valued  functions  s.t.  for  every  f  e  D,  Lf  is  again 
a  real-valued  function,  and  if  f  _>  0  throughout  its  domain, 
then  Lf  >_  0.  Then  L  is  a  positive  linear  operator. 

Remark;  Both  notations  Lf(x)  and  L(f;x)  will  be  used  to  denote  the 
function  Lf  evaluated  at  x. 

Definition  2.2:  If  f  e  C[a,b],  the  modulus  of  continuity  of  f 

on  [a,b]  is  defined  by 

o)f(6)  =  sup{  |  f  (t)-f  (x)  |  |x,t  e  [a,b],  [  x— 1 1  £  5}  for  6  >  0 

Remark;  w  is  nondecreasing  on  (0,°°),  lim  wf(S)  =  0,  and 
f  <5  +  0+ 

WfCyd)  <_  (y+l)w^(6)  for  all  y,6  >  0. 

The  operators  to  be  considered  will  be  of  a  certain  type 
which  R.  DeVore  [2]  calls  interpolation  operators. 
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Definition  2.3:  Let  L  be  a  linear  operator  from  C[a,b]  to  C[a',b’] 

s.t.  there  exist  x  , • • • ,x  e  [a,b]  and  g_,,#,,g  e  C[a',b']  so  that 

in  in 

for  every  f  e  C[a,b]  and  t  e  [a’,bf], 

n 

Lf(t)  =  l  gk(t)f(xk). 
k=l 

Then  L  is  called  an  interpolation  operator,  or  a  proper  interpolation 
operator  if  [a',b']  c  [a,b].  Also,  x^,***,x  will  be  called  the 
nodes  of  L. 

Essentially,  then,  an  interpolation  operator  is  one  for  which 

the  image  of  a  function  is  determined  by  the  values  of  the  function  on 

a  certain  fixed  finite  set  of  points.  Examples  of  interpolation 

operators  are  provided  by  the  processes  of  Lagrange,  Hermite-Fejer, 

extended  Hermite-Fejer,  and  other  forms  of  interpolation.  It  should  be 

noted  that,  if  L  is  a  proper  interpolation  operator,  it  will  not 

necessarily  interpolate  the  values  of  every  function  at  the  points 

x,  ,***,x  (this  occurs  only  when  g,  (x.)  =  6  .  for  all  k,j); 

I  n  k  2  kj 

however,  most  of  the  operators  considered  here  will  do  so. 

The  following  theorem,  due  to  B.  Mond  and  0.  Shisha,  will 
be  used  in  determining  the  rates  of  convergence  for  certain  interpolation 
operators  in  later  chapters  and  represents,  in  a  sense,  the  center  of 
this  thesis. 

Theorem  2.1  (0.  Shisha  and  B.  Mond  [9]):  Let  L  be  a  positive 

linear  operator  from  C[a,b]  to  C[c,d]  whose  domain  includes  the 

2 

functions  1,  t,  and  t  ,  with  LI  =  1.  If  [c,d]  c  [a,b]  and 
||*||  denotes  the  sup  norm  on  [c,d],  then,  for  every  f  e  D, 


. 


5. 


1 1  I*f — f  1 1  co  1  2wf(y),  where  y  =  ||  L  ((x-t) 2  ;x)||  and 

0)^  is  the  modulus  of  continuity  of  f  on  [a,b]. 


2  2 

Remark:  In  forming  L((x-t)  ;x)  ,  L  is  applied  to  (x-t)  as  a  function 

of  t,  and  the  image  is  evaluated  at  x. 


Proof  of  Theorem  2.1.  Let  x  e  [c,d],  t  e  [a,b],  and  6  >  0.  Then 

|  f  ( t )  — f  (x)  |  <_  (l  +  (x-t)^6  (trivial  if  |x-t|  £  6,  and  if 

| x-t  |  >  6  then  |f(t)-f(x)|  _<  ( | x— 1 1  6  ^6)  <  (l  +  |  x— t  [  6  ^)w^(6)  <_ 

2  —  2 

(i  +  |  x-t  |  6  )  0)^(6).)  Applying  L  to  both  sides  of  this  inequality 

(holding  x  fixed)  and  evaluating  the  results  at  x  yields: 

|Lf(x)-f(x)|  _<  (l  +  L((x-t)^;x)5  ^)m^(6). 

The  theorem  then  is  proven  by  putting  6  =  y. 

A  restatement  of  this  theorem,  in  the  context  of  interpolation 
operators,  follows: 


Corollary  2.1:  Let  x, , • • • ,x  e  [a,b],  and  let  m. , • • • ,m  be  nonnegative 
- * -  In  In 

continuous  functions  defined  on  a  subinterval  [c,d]  of  [a,b]. 


n 


I  mk(x)(x-xk) 

Let  p  =  sup  k  -  ,  and  define,  for  f  e  C[a,b], 


xe[c,d] 


n 


I  Mx) 

k=l  k 


n 


l  dl  (x)f(x.) 
k=l 


n 


l 

k=l 


Lf(x)  = 


. 


♦  ts- 


6. 


L  is  then  a  positive  linear  operator  and  if  f  £  C[a,b]  with  modulus 
of  continuity  and  ||  •  ||  denotes  the  sup  norm  on  [c,d],  then 

llLf  -  f  II  oo  1  2wf(v^). 


Remark :  Notice  that  L,  defined  in  this  way,  is  a  proper  interpolation 

n  m^(x) 

operator  -  that  is,  Lf (x)  =  l  g^(x)f(x^),  with  g^(x)  =  -  • 

k”1  I  •.  (X) 

j-1  3 


Moreover,  Lf (x^)  =  f (x^)  whenever  rn^(x^)  =  5  for  all  k,j. 

A  simple  extension  of  this  corollary  will  be  used  to  extend  a 
result  in  Chapter  IV. 


Theorem  2.2:  Let  x, ,***,x  e  [a,b]  and  let  m. , • • • ,m  be  nonnegative 
-  in  in 

continuous  functions  defined  on  [c,d]  <=  [a,b].  Let  ||  •  ||  denote  the 
sup  norm  on  [c,d]  and  let  y  >  1.  Let  L:C[a,b]  ->  C[c,d]  be  defined 


by 


n 


Lf(x)  = 


I  mk(x)f(x) 
k=l  K  * 


n 


I  “1.00 

k=l  k 


Then 


Lf-f  ||  <  2uc 

oo  —  f 


II L  ( I  x-t  |  Y;x) 


Hr 


The  proof  of  this  theorem  is  essentially  similar  to  that  of  Theorem 


2.1. 


■ 
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CHAPTER  3 


HERMITE-FE JER  AND  RELATED  INTERPOLATION  PROCESSES 


In  this  chapter  various  interpolation  processes  and  their 
basic  polynomials  will  be  introduced.  This  information  will  be  used 
in  Chapter  5  to  construct  certain  interpolation  operators  which  will  be 
shown  to  yield  a  good  degree  of  approximation. 

Definition  3.1.  A  polynomial  interpolation  process  on  an  interval 

[a,b]  is  defined  to  be  a  sequence  {L^}  of  interpolation  operators 

such  that  for  every  n  =  1,2, •••  there  exist  polynomials  g  1,»**,g 

nl  nn 

and  nodes  x  <  •  •  •  <  x  _  so  that  for  f  e  C[a,b], 
nn  nl 

V(X)  =  J,  8nk(x)f(xnk)>  X£ta>b]' 

k=l 

The  g  ^  are  called  the  basic  polynomials  of  the  process. 

A  number  of  examples  of  polynomial  interpolation  process  are 
listed  below,  and  the  basic  polynomials  are  given  in  each  case. 

Lagrange  Interpolation:  Lagrange  interpolation  is  perhaps  the  simplest 

form  of  polynomial  interpolation  by  which,  for  an  integer  n  _>  1  and 

a  continuous  function  f (x) ,  a  polynomial  pen.  is  determined  s.t. 

n-1 

p(x  ,  )  =  f (x  .  )  where  x  ,  ,***,x  are  the  distinct  fixed  nodes.  That 

nk  nk  nl  nn 

is,  one  defines  the  Lagrange  polynomial  interpolation  process  {L^} 

with  nodes  x  -,•••, x  in  [-1,1]  so  that  if  f  e  C[-l,l],  then 

ni  nn 

L  f(x.)  =  f(x,),  k  =  l,***,n  and  L  f  e  tt  -  .  The  basic  polynomials 
n  nk  nk  n  n-i 

of  this  process  are 

-  7  - 


■  ■  I  I  ,1  :  I 


8. 


w  (x) 

5,  (x)  =  - - - 

nx  (x-x  .  )w’ (x  ) 


n 


,  where  <o  (x)  =  II  (x-x  .  ) , 
nk'"n'"nk'  n  k-1  nk 


since  then  £  ,  (x  .)  =  6,  .  and  £  .  e  it 

nk  nj  kj  nk  n-1 


Hermite-Fej er  (H-F)  Interpolation:  The  H-F  process  {L^}  with  nodes 


(x  } ,  -1  <  x  <  • • • 
nk  —  nn 


<  x  -  <  1  (n  =  1,2,»**),  is  defined  by: 
ni  — 


n 


Lnf(x)  =  l  k^^(x) f (xnk^  for  f  e  C[-l,l],  x  e  [-1,1], 


k=l 


where  h  ,  (x)  = 
nk 


o>"(x  ,  ) 

1  -  (x-x  )  —77 - 7 

nk  co  (x  .  ) 
n  nk 


?  w  (x) 

£  (x)  ,  £  (x)  =  7 - 7-77 - 7  , 

nk  nk  (x-x  .  )ar  (x  .  ) 

nk  n  nk 


n 


and  co  (x)  =  II  (x-x  .  )  ,  so  that  L  f  e  it. 

n  ,  ,  nk  n  2n-l 

k=l 


Simple  calculations  reveal  that,  for  all  k,j, 


h  (x  .)  =  6 
nk  nj  kj 


L  f  (x  )  =  f(x  ) 
n  nk  nk 


h’  (x  )  =  0  =>  (L  f)  ’  (x  )  =  0. 
nk  nj  n  nk 


Extended  H-F  Interpolation:  This  process  extends  the  previous  one  by 
requiring  also  that  the  second  and  third  derivatives  of  L^f  vanish  at 
the  nodes  xnk>  and  as  defined  by: 


n 


n 


L„f(x)  =  l  Ank(x)f  (x^). 


n 


k=l 


nk" 


where 


Xnk(x)  =  (l  +  7  (x-xnk)  +  C2(x-xnk)Z  +  C3(x-xnk)3)^k(x), 


C  =  -4£’  (x  ) 
1  nk  nk 


. 


9. 


C0  =  10  (x  ,  ))2  -  21",  (x  ,  ) 
2  v  nk  nk  J  nk  nk 


=  _20^nk(xnk^3  +  10llt(xnl)lnl(x.v)  -  kikv> 


nk  nk  nk  nky  3  nk  nk" 


Once  again  calculations  yield  that 


L  f  6  it 
n 


4n-l 


5 


L  f (x  .  )  =  f  (x  )  and 
n  nk  nk 


(V)’(xnk)  "  (Lnf>"(x„k>  “  (Lnf>'"  (xnk>  =  °‘ 


Remarks : 


1) 


is  a  trivial  con¬ 


sequence,  in  each  case,  of  the  properties  imposed  on  the  basic  polynomials 
and  their  degrees  as  polynomials. 

2)  It  will  be  the  task  of  Chapter  5  to  use  information  that 
has  been  discovered  relating  to  certain  special  cases  of  polynomial 
interpolation  processes  in  order  to  derive  "good"  methods  of  linear 
rational  approximation.  There  follows  a  description  (in  some  generality) 
of  the  underlying  methods  to  be  used  in  each  case. 


The  problem  is  to  find  certain  polynomials  g  (x) ,  and 

nk 

n°deS  xnk  (_1  —  Xnn  <  *  *  *  <  xni  —  1} »  n  =  1,2,--*,  1  <  k  <  n,  s.t., 
for  some  interval  I  c  [-1,1], 


(i)  Sn^(x)  il  0  for  all  n,k  and  x  e  I 
(ii)  there  exists  c  >  0  so  that 


sup 
x  e  I 


n 

c  <  )  g  .  (x)  for  all  n  and 

~  k-i  nk 


n 


k=l 


gnk(x)(x-xnk)  I  =  o(1) 


X  €  I 


and 


(iii) 


as  n  ->  °°. 


.* 


. 


' 


10. 


A  sequence  {L  }  of  linear  operators  on  C[-l,l]  whose 
images  are  rational  functions  will  be  defined  by 


n 


l  8nk<*>f 


L  f(x)  = 
n 


k=l 


nk" 


l  8nk<x) 

k=l 


The  Korovkin  method  (Corollary  2.1)  will  be  then  applied  to 
these  operators  to  get  an  estimate  for  ||  L  f*-f  ||  ^  . 


CHAPTER  4 


A  GENERAL  RATIONAL  INTERPOLATION  APPROXIMATION  RESULT 


The  main  results  of  this  chapter  establish  certain  properties 

of  the  functions  h  (x)  (not  necessarily  polynomials),  n  =  1,2,*#*, 

nk 

k  =  l,»**,n,  which  ensure  that  certain  approximation  procedures 
involving  these  h  ^  yield  reasonably  good  degrees  of  approximation. 

Definition  4.1:  Let  [a,b]  be  a  fixed  interval  and  suppose,  for 

n  =  1,2,* ••,  that  points  x  .  (k  =  l,***,n)  are  given  and  there 

nk 

exist  n  >  0  and  integers  nQ,  N  s.t.  for  every  interval  I  c  [a,b] 

of  length  £ (I)  and  for  every  n  >  n  , 

—  o 


card  (klx  .  e  1}  <  n  •  £(I)  •  n  +  N; 
1  nk  — 


then  the  points  {x  are  said  to  be  Mn-dispersed  in  [a,b]" 


Theorem  4.1:  Let  points  -1  <  x  <  •••  <  x  .  <  1  and  functions 
-  —  nn  nl  — 


h  - (x) , *  *  * ,h  (x)  be  given  (n  =  1,2,***),  and  suppose  there  exist 
nl  nn 

constants  n  >  0,  M>0,  e  >  p  >  0,  c2  j>  C.^  >  0,  and  an  integer  ^ 

s.t.  for  all  x  e  I  =  [-l+e,l-e],  the  following  conditions  hold: 

e 

(1)  h  ^(x)  >_  0,  all  n  >_  n^,  all  k  =  l,***,n 


n 


(2)  ci  1.  1  hnk(x)  —  C2’  a11  n  21 

k=l 


(3)  h  ,  (x) 


M 


nk  —  2 ,  v 

n  (x-xnk) 


—  ,  all  n  _>  n^,  all  k  =  1, 


,n 


11 


: 


■ 


12. 


(4)  The  points  {x^}  are  ^-dispersed  in  1^  =  [-l+y,l-y]. 
Define  a  sequence  of  positive  linear  operators  on  C[-l,l]  by: 


n 


Lnf(x)  = 


kIihnk(x)f(xnk) 

n 


x  e  I 


2  hnk(x) 


k=l 


Then,  if  II  *  IL  denotes  the  sup  norm  on  1^, 


IIVfll»  -c“i 


A.n  £n  n 


n 


where  C  is  a  constant  independent  of  n  and  f,  and  denotes  the 

modulus  of  continuity  of  f  on  [-1,1]. 


Remarks :  This  result  is,  in  a  certain  sense,  a  generalization  of 

Theorem  5.1.  The  degree  of  approximation  obtained  in  the  theorems 
of  Chapter  5  is  o(u)f(^-)),  while  the  result  in  Theorem  4.1  is  the 
degree  0 


_ f  n' 

oof  (--^n  ^  -n)  .  The  factor  / £n  £n  n  converges  to  00 


as  n  ->  °°;  however,  numerically,  it  grows  rather  slowly.  For 

1010  7 - 

example,  if  n  =  10  ,  still  /£n  £n  n  <  5. 


The  proof  of  Theorem  4.1  requires  some  lemmas: 


Lemma  4.1:  Let  points  x  functions  h  ^  (n  =  1,2,««»;  k  =  l,*»*,n), 
constants  n  >  0,  M>0,  e  >  y  >  0,  and  >  0,  and  an  integer 


n^  be  given,  so  that  assumptions  (1)  to  (4)  of  Theorem  4.1  are  satisfied. 


(i) 


if  a  >  0 
n  _>  n2»  x 


there  exists  an  n ^  and 


A  >0  s.t. 
a 


Then 


' 


. 


13. 


and  (ii)  if  {a  }  is  a  sequence  of  positive  reals  s.t.  a  4-0 

n  n 

and  n  ^  ^  I  £na  I  -y  0,  then  for  some  n_  and  A  >  0, 
1  n1  3  * 


a 


n 


n  >  n„,  x  e  I  =>  n 
—  j  e 


n 


k=l 


£na  |  T  h^.  (x)  >  A. 
n „  nk  — 


Lemma  4.2:  Under  the  assumptions  of  Lemma  4.1,  there  exists  a  constant 

B  >  0  and  an  integer  n  s.t.  whenever  n  >  n, ,  x  e  I  ,  we  have 

4  —  4  e 


l  hnk(x)(x_Xnk)2 

k=l  n  _  <  B  £n  £n  n 


n 


l  hnk(x) 

k=l  nK 


n 


Remarks :  Part  (i)  of  Lemma  4.1  is  actually  implied  by  part  (ii) ,  but 

the  statement  and  proof  of  part  (i)  are  given  since  they  form  the 
basis  for  part  (ii) ,  and  the  proof  is  simpler  in  the  former  case. 


It  should  be  noted  that  essentially  the  best  choice  for  the 

a  is  a  =  (£n  n)  -  i.e.,  this  choice  gives  essentially  the  slowest 

a 

growth  in  the  factor  n  Una  I.  Suppose  (g  }  is  any  other  choice 

nn  n 

g  a 

s.t.  g  >0,  g  -*  0,  and  n  n|£ng  I  <  Kn  n  I  £na  I  for  some  K  >  0; 
n  n  0  '  n 1  —  1  n 1 

g  a 

then  I  £ng  I  -*  00  =>  n  n  =  0(£n  in  n)  (since  n  n  I  £na  I  =  £n  £n  n)  ==> 
n  n 


g  =  o(- 

n  v 


£n  £n  £n  n^ 


£n  n 


£ng  >  —  £n  £n  n  for  n  large  enough,  and 
n  —  L 


g. 


hence  n  n|£ng  I  >  £n  £n  n  for  n  large  enough.  Therefore  (except 
n  —  L 

for  at  most  a  constant  factor),  a  =  (£n  n)  ^  yields  the  best  result. 

n 


From  the  above  it  is  also  clear  that  the  requirement 

n  ^  |  £na  I  o  does  not  weaken  the  result  of  Lemma  4.1,  since 

n 

certainly  n  |  £na  |  -y  q  for  the  choice  a^  =  (£n  n) 


Proof  of  Lemma  4.1:  It  can  be  assumed  w.l.g.  that  j>  1.  (If  not. 


' 


■ 
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+  -1+2 

divide  by  C^.).  Let  n  e  Z  and  define  v^=n  ,  i  =  1,2,*** 

2  + 
so  that  nv.  =  v.  and  v.  ,  <  v.  for  all  i.  Let  n  ,  N  e  Z 
i  l-l  l-l  l  o y 

be  as  in  Definition  4.1  (on  the  interval  [a,b]  =  [-1+y ,1-y ] . ) . 


Part  (i) :  Let  ot  >  0,  and  assume  w.l.g. 


j  e  Z  s 


. t.  2  ^  <  a  <  2  j+1. 


Let  n~  >  n  s.t.  n  >  n~ 
2—o  —  2 


Let  x  e  I  and 

Cn-d/2)  <  e  _  y# 


I1(x)  =  (k 


x-x  .  >  Cv  } 

1  nk 1  1 


I  GO  =  {k 

m 


Cv  <  x-x  ,  <  Cv  .. 

m  '  nk 1  —  m-1 


!j+1(x)  =  {k 


x-x  .  <  Cv  .  } 

nk1  -  j 


that  a  <  1.  Let 
C  =  max(l,  2j(2n+N)M}. 
Assume  n  >_  n^.  Define 


},  m  =  2 ,  •  •  • ,  j 


Then,  for  m  =  2,***,j+l, 

o  o 

card  fl  (x)l  <  n  •  n  •  2Cv  +N  <  (2Cn+N)nv  ..  =  (2Cn+N)n  v 
v  m  J  —  m-1  —  m-1  m 

By  the  above  and  assumption  (3)  ,  for  1  <_  m  _<  j  ,  simple  calculations 
reveal  that 


m 


so  that 


l  hnk«  1  I  • 


ked+1(x) 


and  hence 


-2 


"J 


^  ^hnk(x)  -  4  card[lj+1(x))  -  8Cn+N  -  A«n 


kel^.+1(x) 


with 


A  = 


a  8Cn+4N 


' 


. 

■ 
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Part  (ii) :  The  proof  of  part  (ii)  is  much  the  same  as  that  of 

.  _+ 


part 

(i) .  For  each  n 

s.t. 

a  <1,  choose 
n 

=  0(|  inot  | )  as 
n 

2A 

Jn+1 

<  a  <2  .  Then  j 

—  n  Jn 

for 

n  large  enough) . 

Also 

let  C  =  max(l, 
n 

C  < 

n  — 

K  Jin  a  for  some 

n 

o 

A 

independent  of 

n 


n 


n 


—  o 

n  >  m  =>  C  n  <  e  _  y  (possible  since  n  ^^)  Una  I  ->  0) . 

—  n  r  1  n1 

Now,  as  in  the  proof  of  (i) ,  if  n  >  m  and  x  e  I  ,  then 

—  e 


^  hnk(x)  -  ^8C  n  +  4N^n 
k=l  n 


Let  n^  2.  m  s.t.  n  n^ 


a  <1  and  C  n  >  4N  (possible  since 
n  n 

n  2 

a  0  =>  j  ->  <»  =>  c  ->  °°)  and  then  n  >  n„  =^>  j  h  .  (x)  > 
n  n  n  j  .  1  nR 

k=l 

-a 

n  n  1  an  -1 

~9C — n  —  9Kn  l£nanl^  *  Hence  (ii)  is  established  with 
n 


A  = 


9Kn  * 


Proof  of  Lemma  4.2:  This  follows  immediately  from  Lemma  4.1  by  setting 
-1 


=  (Jin  n)  ,  which  yields,  for  n  >_  n^. 


B  2  2  an  i  |  -2  p 

I  h  .  (x)  (x-x  .  )  n  £na  Mn  J  h  ,  (x) 
nk  nk  1  n 1  ,  nk 

k=l _  _ k=l _ 


n 


2  hnk(x) 


k=l 


(Jin  Jin  n)  •  MC, 

- A  " 


MC, 


so  the  lemma  is  established  with  B  = 


— —  and  n.  =  n_. 
A  4  3 


Proof  of  Theorem  4.1:  By  lemma  4.2,  we  have 

2  n  _  B  Jin  Jin  n 


Ln((x-t)  ;x)  £ 


n 


: 


. 


, 
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Theorem  4.1  follows  immediately  from  this  inequality  and  Corollary  2.1. 

A  few  easy  modifications  in  the  proofs  of  the  above  lemmas 
verify  the  following  extension  (see  Theorem  2.2  and  the  following  remark). 


Lemma  4 

B  >0 
P 


.  3 :  If  p  is  an  integer,  p  _>  2,  then  there  exists  a  constant 

and  an  integer  m  s.t.  whenever  .  n  >  m  ,  x  e  I  ,  we  have 
&  p  —  p  e 


\hnk(x)  (x-Xnk)2(P~1) 
k=l _ 

l  hnk« 

k=l  nK 


<  B 


fa  fe  n 


n 


P-1 


Using  this  lemma,  one  can  easily  prove  a  theorem  which  actually  implies 
Theorem  4.1  with  the  choice  p  =  2: 


Theorem  4.2.  Let  p  be  any  integer  with  p  _>  2,  and  let  points 

-1  <  x  <***<x1<l  and  functions  h  -(x),***,h  (x)  be  given 

—  nn  ni  —  ni  nn 

for  every  n  =  1,2,***;  suppose  there  exist  constants  n  >  0,  M  >  0, 
e  >  y  >  0,  C2>C^>0,  and  an  integer  n^  s.t.,  for  x  €  I£, 
assumptions  (1)  thru  (4)  in  the  statement  of  Theorem  4.1  are  satisfied. 
Define 


for  f  e  C[-l,l],  x  e  1^. 

Then,  if  || • denotes  the  sup  norm  on  I£,  there  exists 
0^  >  0  independent  of  n  an<3  f  s.t. 

/ £n  5,n  n' 


||l  f-f  ||  <  c  w_ 

11  n.p  " 00  —  p  f 


' 


. 
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The  following  result  will  assist  in  demonstrating  the 
applicability  of  Theorems  4.1  and  4.2  to  several  examples: 


Theorem  4.3.  Let  p^(x)  be  one  °f  the  classical  polynomials  T  (x) , 
Pn(x)>  ^(x),  ^n+l^X)’  °r  ^-x  )P^_-^(X)>  an^  ^et  Pn  have  zeros 


-1  <  x  <  •••  <  x  -  <  1,  for  n  =  1,2, 
—  nn  nl  —  ’ 


Then  for  every  y  >  0 


there  exists  n  >  0  s.t.  the  x^k  are  n-dispersed  in  [-l+y,l-y]. 

Proof:  It  is  clearly  enough  to  prove  this  for  the  choices  p  =  T  or 

-  n  n 

P  ,  since  the  zeros  of  the  other  choices  alternate  with  the  Tchebycheff 
or  Legendre  zeros,  by  the  Mean  Value  Theorem.  T^  and  P^  are  examples 

( C X  C l) 

of  the  classical  ultraspherical  polynomials  P^  ’  (x) ,  with  a  =  -  y  , 

0  respectively.  Therefore  a  separation  theorem  (Szego  [10],  pp.  122) 

due  to  Stieltjes  can  be  applied,  and  states  that  if  x^^  =  cos  0^k, 

,  .  ,  _  „ (a,  a)  .  .  1  1 

k  =  l,***,n,  are  the  zeros  of  P  (x)  with  -  —  <  a  <  — ,  then 

n  Z  —  —  Z 


(k  -  <  0  v  <  k 

^  2  n  —  «+  — 


nk-"nll’  k  =  <!.«..  6nk  <  0) ; 


for  the  k  for  which  0^k  >0,  a  similar  result  clearly  occurs  by  the 

IT 

symmetry  of  the  zeros  about  the  point  0  =  —  • 

Let  y  >  0;  then  there  exists  some  y  >  0  s.t.  x,t  e  [-l+y,l-y] 


.n. 


=>  I  x-t  I  _>  y  I  arc  cos  x  -  arc  cos  1 1  .  .  *  .  if  1  k  _<  [— ]  -  1  and 

x  .  ,  x  ,  -  e  [-l+y,l-y],  then 
nk  nk+1 


c  ,  -  x  .  >  y  arc  cos  x  ,  -  arc  cos  x  , 

nk  nk+1 1  —  1  nk  nk+1 


Y'9nk  9nk+l 


r(k  +  t) 


>  Y 


kTT 


n 


n+1 


YTT 

2n  * 


and  similarly  for  pairs  x  k>  xnk+1  on  the  other  side  of  zero 

(k  =  [~~2~ ]  >  *  *  *  »n)  •  Two  cases  are  considered  (n  even  and  n  odd)  for 

the  zeros  at  or  next  to  the  point  zero.  If  n  is  even,  let  k  =  ^  =  [-j]  , 


' 


■ 


■ 
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and  then 


„  kir  (k+1)  tt 
<  — —  <  — - - —  < 


6nk  —  n+T  iSf  -  9nk+l  (by  symmeCry) 


Xnk  Xnk+1  ’  Y 


yrr 

2n  * 


If  n  is  odd,  let  k  =  =  [y] ,  and  then 


e  <  Jil_  <  1  =  e  <  tt+jgjhr  <  e 

nk  —  n+1  2  nk+1  n+1  —  nk+2  2  n+1 


TT  _  (k+1)  IT 


so  that  x  ,  -x  ,  , , 

1  nk  nk+1 


=  x  ,  . , -x 


hk+1  nk+2 1  —  n+1  2n 


IT  YTT 

>  y  — —  >  J 


Therefore  the  distance  between  any  two  consecutive  zeros  x  ^  and 

YTT 

x  .  in  [ — 1+y , 1— y ]  is  at  least  tt-  •  Therefore  if  I  is  any 
nk+1  2n 

subinterval  of  [ — 1+u , 1— y ]  of  length  £(I)  then 


card  {k 


X  .  €  1}  <  —  •  £(I)  +  1, 

nk  —  yir 


so  the  {x  ,  }  are  p-dispersed  in  [-l+y,l-p]  where  n  =  —  . 

nk  Y77 

( ot  ct) 

Remark .  The  classical  ultraspherical  polynomials  *  (x)  satisfy 

a  differential  equation  of  the  form 


(l-x^)y"  -  vxy'  +  n(n  +  v  -  l)y  =  0, 


where  v  =  2a  +  2.  Note  that  if  p  satisfies  such  an  equation  and 

n 

2 

has  roots  x  ,  (k  =  l,***,n),  and  if  q  (x)  =  (1-x  )p  (x) ,  then 
nk.  n  n 


Pn" (xnk>  V  Xnk 

P  *(x  ,  )  .  2  ’ 

n  nk  1-x  , 
nk 


q"  (x  ,  )  (v-4)x 

n  nk  _  _ nk 

q’ (x  )  2 

n  nk  1-x  . 

nk 


and 


■ 


•  . 

. 


Therefore 


(i  - 


p"(x  .  ) 

(x-x  )  ■■  ~  -  -n  ■-) 
nk  Pn^nk5 


1  -  VX  ,  X  +  (v-l)x  , 
nk  nk 


1  -  x 


nk 


and 


(i  - 


q"(x  .  ) 

,  v  n  nk  'v 

(x-x  )  — - v-j 

nk  q  (x  ) J 
n  nk 


1  -  (v-4)x  x  +  (v-5)x 

nk  nk 


1  -  x 


nk 


In  each  of  the  following  examples,  -l<x  <  •••  <  x  ,  <  1 

—  nn  nl  — 

will  denote  the  zeros  of  p  (x) .  It  should  be  noted  that  Theorem  4.3 

n 

serves  to  establish  in  each  case  that  the  x  ^  are  p-dispersed  in 
[-l+y,l-y]  for  each  y  >  0  and  appropriate  p.  It  remains  therefore 
to  verify  assumptions  (1) — (3 )  of  Theorem  4.1  in  each  case. 


Example  1; 


P_(x)  =  T  (x),  h  (x)  = 
n  n  nk. 


(1-x  x)T  (x) 
nk  n 


(1-x  x)T  (x) 
nk  n 


(1"Xnk)(x"X-^)2^(x-^^2  n2(x"X^)2 


nk  k  n  nk 


nk 


Hence 


0  <  h  .  (x)  < 


nk  —  2  (  s  2 

n  (x-xnk) 


n 


and  \  h  (x)  =  1  (see  the  remarks  at  the  end  of  Chapter  3,  as  well 
k=l  nk 

as  the  definition  of  H-F  interpolation.)  In  this  case,  one  can  choose 


any  e  >  y  >  0. 


Example  2;  Pn(x)  =  (x) »  (x)  = 


n 


nk 


a  -  2xnkx  +  x^k)P^(x)  ^ 

a-x2j(x-xnk)2(p;(xnk))2 


nk 

n 


Once  again  (x)  ^0  on  [-1,1]  and  \  hnk^  "  1*  The  ProPerty 


nk 

that  for  e  >  0  there  exists  an  M  s.t. 


k=l 


hnk(x)  -  “2 


M 


n  (x-xnk) 


■ 


. 
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for  all  x  £  I  ,  and  all  n,  k,  requires  certain  estimates: 


Pn(x)|  =  0(n  (1/2)(l-x2)  (1/4)),  all  n 


1  "  Xnk  " 


,2  “2 

k  n 


n 


,  if  k  _<  —  ,  all  n,k 


(n+l-k)2n  2,  if  k  >  y 


ak  (3/2)n2  <  |p’(x  )|  <  a  k  (3/2)n2,  if  k<^ 

a  (n+l-k)"(3/2)n2  <  |p' (x  ,  ) |  <  a  (n+l-k)" (3/2)n2 ,  if  k  , 

A  —  n  nk  —  i  A 

with  constants  a^  >  a^  >  0.  These  estimates  may  be  found  in  Szego  [10; 
pp.  122,  165,  238], 

Let  e  >  0  and  restrict  x  e  1^.  Denoting  j  =  min{k, n+l-k} , 
it  follows  from  the  above  inequalities  that  there  exists  K  >  0 
independent  of  n,k,x  s.t. 


hnk(x)  - 


4K 


. 2  -2  r.  2  .-3  4 ,  .2 

2  n  n^l-x  j  n  (x-x_.,  ) 


nk 


=  0 


n2(x-xnk)2 


Example  3:  p  (x)  =  U  (x) ,  h  .  (x)  = 
- -  ’■*  n  nk 


n 


n 


(l-3x  .x  +  2x2.  )U2(x) 
_ nk _ nk  n _ 

(1'Xnk)(x“X">')2(U"(X"t)^2 


nk  ^  n  nk' 


3-2/2 


In  this  case  it  is  necessary  to  choose  e  _>  — 3 —  to  ensure  that 

n 


x 


£  I  =>  h  ,  (x)  >0.  We  do  have  T  h  .  (x)  =  1. 
c  nk  — 


k=l 


nk 


|u’(x  .)|  =  and  U2  (x)  <  sin  2  (arc  cos  x)  =  0(1) 

n  nk  ,2  n  — 


1-x 


nk 


Since 


' 

■ 
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on  I  ,  it  follows  that 
e* 


hnk«  ^ 


6(l-xnk>  •  0(1) 

(n+l)2(x-xnk)2 


=  0 


'I'2(x-Xnk) 


<i-Ax„vx  +  3x!i,>  w)2 


Example  4:  Pn(x)  =  Pn+i^x^»  hnk(x)  =  2 


nk 


nk  v  n+1 


Cl-Xnk)(x-Xnk)2fPn+l(Xnk)52 

2-/3 

If  e  =  - ^ —  ,  then  x  e  I  =>  h  .  (x)  >  0  for  all  n,k. 

z  e  nK  — 

n 

Also,  again,  £  h  ,  (x)  =  1. 


k=l 


nk 


Two  inequalities  to  be  used  can  be  found  in  Balazs  and  Turan  [2]  or 
Szego  [10] : 

|P’n+1(x)|  1  /2  (n+2)  (l-x2)"(3/4) 


and 


|Pn+l(xnk)l  -  (8,,j) 


-(1/2) 


where  j  =  min {k, n+2 -k}. 

Furthermore,  the  differential  equation 

(l-x2)P”+1(x)  -  2xP^+1(x)  +  (n+l)(n+2)Pn+1(x)  =  0 

yields 


(1-Xnk)Pn+l(xnk)  =  -(n+1) (n+2)Pn+l(xnk) ' 


so  that 


hnk(x)  - 


8(l-x2  )  •  2 (n+2) (1-x2)  (3/2) 
_ nk _ 

(x-x^)  2  (n+1)  2  (n+2)  2 


=  0 


2  f  .2 

n  (x-xnk) 


•  8tTJ 


. 


' 


22 


?  (1-x2,  )(d-x2)p;+1  (X))2 

Example  5;  Pn(x)  =  (l"x  )p,!1+1(x)>  (x>  =  - J~ 


n+1 


(1'Xnk)  (X-Xnk)2^(xnk^2 


In  this  case  the  {x  ,  }  are  taken  to  be  the  zeros  of  P',(x)  (the 

nk  n+1 

other  roots,  x  =  ±1,  of  p  (x)  are  disregarded  for  simplicity). 

n 


Letting 


h  (x)  = 
no 


p"(l) 


n 


1  -  (X-1)  -y 


p;a>j 


2,  \ 

pn(x) 


(x-l)2(p'(l))2 

n 


and 


h  ,  (x)  = 

n,n+l 


p"  (-l)'i 


n 


1  -  (x+1) 


Pn(K) 


(x+l)2(p'(p'(-l))2 
n  n 


one  gets 


n+1 

E  hnk(x) 


k=o 


=  1. 


Now,  p’ (x)  =  -(n  +n)P  (x)  (from  the  differential  equation) 
n  n 

2 

and  p"(x)  =  -(n  +n)P’(x),  so  that 
n  n 


Pn(±1)  I  =  i(n2+n)» 


p"(±l) 


n 


P;<±i) 


n  +n 
2  * 


Therefore  if  k  =  0  or  n+1,  £  >  0,  and  x  e  I  , 

e 


lhnk(x)l  - 


1  + 


2(n2+n)l  (U-x2)P;(x))2 


e2 (n2+n) 2 


=  0(4)  ((l-x2)p;(x))2 


n 


=  0  f— 1  ,  using  |P*  (x)  |  <  /2  (n+1)  (1-x2) 

1  n  1 


Therefore,  for  n  large  enough,  |hnQ(x)  +  I  <  2  for  a11 


x  e  I  ,  so  that 
e 


- 


' 
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1 

2 


< 


l 


k=l 


hnk(x) 


Finally,  if  1  _<  k  n  and  x  e  1^,  then 

(1‘Xnk)(1'x2)2(Pn(x)^2 

h  ,(x)  =  - r - r - - - o  0 -  °»  clearly> 

(1-x  )  (x-x  )  (n  +n)  P  (x  ) 

nk  nk  n  nk 


((l-x2)3/4P^(x)]2vCx^ 

2  2  2  2 
(n  +n)  P  (x  ) (x-x  ,  ) 
n  nk  nk 


2(n+l)vl-x  (8-n-j) 

>2 


(n2+n) 2 (x-x 


nk 


where  j  =  min{k,n+l-k}  (using  the  two  inequalities  quoted  in 
Example  4) 


=  0 


2  ,  v  2 

n  (x-xnk) 


From  these  examples  it  is  evident  that  Theorems  4.1  and  4.2  apply  to 


a  large  number  of  cases. 


CHAPTER  5 


ACHIEVING  A  BETTER  DEGREE  OF  APPROXIMATION 


In  this  chapter,  three  approximation  theorems  which  are 
related  to  the  results  of  Chapter  4  will  be  stated  and  proven.  The 
special  properties  of  each  sequence  of  operators  will  be  used  to  attain 
good  degrees  of  approximation. 


As  in  the  previous  chapters,  will  denote  the  modulus  of 

continuity  of  f  e  C[-l,l]. 

Theorem  5.1.  Let  x  , (-1  <x  <***<x.<l)  be  the  zeros 
-  nl  nn  nn  nl 

of  T  (x)  and  define  L  :C[— 1,1]  C [  — 1 , 1  ]  by 
n  n 


n 


l  mnk(x)f(xnk) 


L  f (x)  = 


k=l 


n 


^  mnk(x) 


k=l 


where 


(x)  ?  ?  1  ? 

mnk(x)  =  - /  1 (1’X  > (1'Xnk)  +  2(x'Xnk)  1 ' 

n  (x~xnk) 


Then  L  is  a  positive  linear  interpolation  operator,  and  if 
n 

denotes  the  sup  norm  on  [-1,1],  we  have 


GO 


L  f-f 
n 


<  4w  (— )  . 
—  f  vny 


Proof  of  Theorem  5.1:  In  [6]  there  appears  the  following  representation 
of  the  basic  polynomial  X  ^(x)  of  the  extended  H-F  polynomial  inter¬ 
polation  process  with  Tchebycheff  nodes  (see  Chapter  3) : 


24  - 


■ 


■ 

' 


* 


25. 


Xnk(x)  ’  ~4 


4 

T4(x) 

n 


l(Kix)2  +  (x~x  J 


2  2 (n  -1) 


n  (x"Xnk} 


4  1  “nk“/  v“  “nk/  3 


(1'Xnkx)  -  ¥  XnkX  J 


which  can  be  rewritten: 


\ik(x)  =  “4,  ,4 

n  (x-xnk) 


TJ(x)  ,  ,  ,  ,  .  2 

{(1-x  )(l-xnk)  +2<*-*nk)  +  (*V~ 


)(x-xnk)2a-xnkX)} 


so  that 


n 


1=1  Ank(x) 

k=l 


?  ,  s  ,  (4n2-l)  2.  .3  Tn(x)(1"Xnkx) 

-  J“nk(x>  +  Vx) .  i,  - 72- 

k=l  6n  k=l  n  (x-x  .  ) 

nk 


Tn(x) (l-xnkx) 

Now,  since  - - - —  ,  k  =  l,***,n,  are  the  basic  polynomials 

n2(x-xnk)2 

“  1 

for  the  H-F  process  on  the  nodes  xnk>  one  obtains  2,  m^Cx)  —  T  » 


k=l 


nk 


Also 


n  ?  n  Tn  (x)  9  9  9  9 

!  fflnv(X)Kv)  =  l  “7 - 9  )l 


.  ,  nk  nk  .  4  ,  .2  v  nk 

k=l  k=l  n  (x-x  ,  ) 

nk7 


_1_  ?  Tn(x)(1-xnkx)  J_ 

-  2  1  \  2,  N 2  2  ’ 

n  k=l  n  (x-x  .  )  n 

nk 


and  so  the  result  follows  from  Corollary  2.1. 

Theorems  5.2  and  5.3  deal  with  approximation  by  interpolation 

operators  having  as  nodes  the  zeros  of  U  (x)  and  of  P  (x) . 

n  n 

Theorem  5.2:  Let  x  ,  x  (-1  <x  <***<x,<l)  be  the  zeros 

-  nl  nn  nn  nl 

of  U  (x)  and  let  e  >  0.  Define  L  :C[-1,1]  -*  C(I  )  by: 


n 


n 


n 


L  f(x)  = 
n 


I  n,nk(x)£(xnk) 
k=l 


n 


1  n.nk(x) 


k=l 


y 


i 


' 


26. 


where 


“nk(x)  = 


un(x) 


(x-x  )U' (x  ) 
v  nk  n  nk  ' 


Let  ||.|L  denote  the  sup  norm  on  I£.  Then  there  exists  a  constant 
A  >  0  and  an  integer  n^  s.t.  for  every  f  e  C[-l,l]  and  every 
n  >_  nx, 

|| L  f-f  ||  <  Aw  (-)  . 

"  n  °°  —  f  hi' 

Remark.  The  fact  that  the  error  estimate  in  Theorem  5.2  holds  only  on 
I£,  and  not  on  [-1,1],  as  in  Theorem  5.1,  is  not  really  a  restriction, 

for  if  f  e  C[-l,l]  one  can  approximate  f  on  the  whole  interval 

with  a  modest  amount  of  juggling: 

/ 

f (-1)  ,  x  <  -1  +  e 
Define  g(x)  =  /  f  >  x  e  I 

f(l)  ,  x  >  1  -  e 

V 

so  that  w^(6)  =  w£((l_e)  Ji  ((1-e)  ^  +  l)m^(6).  By  Theorem  5.2, 
if  x  e  Ie, 

|Lng(x)  -  g(x)|  <  Aug^)  <_  A  ((1-e)  1  +  l)wf(^)  , 
and,  letting  t  =  ^7  ,  one  gets  that,  for  t  e  [-1,1], 

| L^g  ( (1-e)  t)  -  f  (t)  |  <_  A  ((1-e)  • 

Thus,  one  obtains,  by  changes  of  variables,  an  approximation  to  f  on 
all  of  [-1,1] . 

The  proof  of  Theorem  5.2  will  require  some  knowledge  of  the  H-F 

and  extended  H-F  polynomial  interpolation  processes  on  the  nodes  x  . 


■ 


' 


27. 


If,  as  in  Chapter  3,  one  denotes  by  h  .  (x)  and  X  ,  (x)  the  basic 

nk  nk 

polynomials  of,  respectively,  the  H-F  and  extended  H-F  process  on  the 
nodes  x^,  then  one  has  the  representations: 


hnk(x) 


(l-3x  x  +  2x  .  ) 
nk  nk 


1  -  x 


Vx> 


nk 


(x-x  .  )U' (x  .  ) 
v  nk  n  nk  ' 


and 


*nk(x)  =  t1  +  Cl(x-xnk)  +  C2(x-xnk)2  +  C3(x’xnk)3) 


u  (x) 
n 


14 


(x-x  )U' (x  ) 
nk  n  nk 


where  (see  Sharma  and  Tzimbalario  [8]) 

-6x 


C,  = 


nk 


1-x 

25x 


= 


2 

nk 

2 

nk 


-  2 


+ 


2N 


C„  = 


2  2 

2(1‘Xnk) 

^1-x2  ' 

nk 

15x  .  -  35x3. 
nk  nk 

ION  x 

3  (1-x 


nk 


2  3 

2(1-Xnk> 


2  2 

3(1-xnk> 


nk 


(N  =  n  +  2n) 


Some  other  facts  will  be  needed: 

(1) 


U'(x  )  = 

n  nk  1 


n+1 


1-x 


nk 


(2)  If  0  <  £  <  1  and  M  =  — : 


e  sin  arc  cos(l-e) 


<  00 ,  then,  for  every 


x  e  I  ,  U  (x)  <  M 

£  n  £ 

<3>  ”  "nk(X)  =  (x-Xnk>U;(xnk> 


for  all  n,k,  and  =  max(M£,2), 


then  there  exists  m  s.t.  if  n>m,  l<k<n,  and  x  e  I  ,  then 

o  —  o  —  —  £ 

111  .  (x)  I  <  K  .  To  see  this  is  not  difficult.  First,  let  m  be  s.t. 

1  nkv  1  —  e  o 


n  >  m 
—  o 


— <  /6  M  n  ,  Let  n,k,x  be  fixed  s.t.  x  e  I  , 

n+1  £  —  12 


n>m,  l<k<n,  and  put  M  =  /6  M  n 
—  o  —  —  £ 


-(1/2) 


x-x 


£  ' 

1 


nk1  —  n+1 


Either 


1  I  I 


28. 


or  x-x  .  <  — —  .  In  the  former  case,  one  has 

nk 1  n+1  ’ 


*'nk(x)l  = 


VX>(1-Xnk> 


(x-xnk) (n+1) 


<  M  . 
--  e 


In  case  x-x  .  <  — — ,  we  have 

‘  nk 1  n+1 


I  .  hni(x)  1  1  +  ,  l  .  Ihni(x> I 

x-x  .  <M  J  x-x  .  >M  J 

nj  1  nj  1 


=  1  + 


x-x  .  >M 
nj 


(l-3x  .x  +  2x2 . ) (1-x2 . )IJ2 (x) 

B2_ BJ  B2  B 


(n+1) 2 (x-x  .  )2 
nj 


<  1  + 


6 • M2 • n 

e 


x-x  .  I  >M  6M2  (n+1) 2 
nj  1 


<  2  . 


Now,  if  | x-x  .  |  <  M  <  -^r  ,  (1-x2  )  >  (1-x  )  >  —  ,  so 

nj  —  —  12  nk  nk  2 


h  .  = 
nj 


1  - 


1  - 


3x  .  (x-x  . ) 
nj  nj 


1-x2. 

i  ,n3 

6  x  .  M'i 
nj 


£2  .  (x) 
nj 


A  ,(x) 

nj 


>  i  £2.(x)  >  0. 

—  2  nj 


\  h  .  (x)  <4  so 


Therefore  £  (x)  <  2h  .  (x)  <  2  • 

nk  —  nk  —  i  “  i  nj 

x-x  .  <M 
1  nj 

£  .  (x)  <  2  <  K  . 
nk  —  e 

Since  the  first  finitely  many  £^k(x)  are  als°  uniformly  bounded  in 


I£* 


K  =  sup { | ^nk(x) I 


n  =  1,2,***,  l<k<n,  xel}<°°. 

—  —  e 


(4)  Let  p  ,  (x) ,  l<k<n,  n  =  1,2,* ••  be  given  s.t. 
nk.  —  — 


B  =  sup{ |pnk(x) | 


l<k<n,  xel}<°°. 
—  —  e 


, 


t|"  „.i  • 5  -  '<£i£3||B  . i 


* 


29. 


Let  n  >  0  s.t.  the  x^  are  n-dispersed  in  I  (e  >  0  fixed). 


Then 


|  Pnk(x)(x-XnkHnk(x) 


k=l  1-x 


nk 


0(£ 


where  the  0  constant  is  independent  of  x  e  I  . 
is  elementary: 


The  proof  of  this 


Pnk(x)l|x-xnk|lnk(x) 


x-x  .  < 

nk 1 


1-x 


nk 


IB  l 


U  (x)  A  .  (x)  M  •  K 

'  n  M  nk  —  <  B  •  0(v£)  •  £ 


n+1 


n+1 


x-x  ,  <■ 

1  nk 1 


and 


I 


Pnk(x) 


X~Xnk^  ^nk(x) 


x-x 


1-x 


nk 


nk 


IB-  l 


x-x 


nk  '  — 


(l-x2,)U2(x)  B 
nk  n  


n  •  M  •  ViT 

e 


1  (n+1)  | x-x 

vk 


nk 1 


(n+1) 


so 


n  (x“*  ? 

l  Pnk(x)  — ^  'nk(x) 
k=l  1-x  , 

nk 


as  claimed. 


„  r„2/N  ,  ,  n/h  .  as  n->°°,  where  the  "0"  constant 

Lemma  5.1:  l  l  ,  (x)  =  1  +  0[— J 

k=l  vC 


is  independent  of  x  e  I  . 


. 

' 


■ 


30. 


Proof : 


Lemma  5 . 2 


uniformly 

Proof : 


where 


n 

so  l 
k=l 


n 


n 


L  =  l  hnk(x) 
k=l  n 


=  l  *nk(X)+  l 


k=l 


n  3x  (x-x  ) 

nk  nk  „2  .  . 

- , -  l  ,  (x) 

i  2  nk 

k=l  1-x 

nk 


n 


£  £^^.(x)  +  0  (-^)  by  (4)  above. 


n 

I 

k=l 


‘  a  2N  U^(x)(1-x2  )£2  (x) 

(x)  +  - 2 - ^L_Elk - 

.  nk  3  (n+1) 2 


for  x  e  I  ,  as  n  ->  00 . 
e 


n 


I  Xnk(x)  =  and 


k=l 


Xnk(x)  =  t1  +  Cl(x-Xnk)  +  C2(x'Xnk)  +  C3(X_Xnk)  ) 


3V 


nk 


-6x 


ci  = 


hk 


1-x 


nk 


2N 


29x  -4 

C2=-JVr  +  -  2 

-  2(l-xnk)  3(l-x  ) 


C3  = 


15x  .  -35x  . 
nk  nk 

2  3 

2(1-Xnk> 


ION  x 


nk 


2  2  * 
3<1-xnk) 


n 


"k(x)  =  Ji 


,  2N (x-x  )  . 

*nk«  +  "  "  t\—  *nk(x) 


3(1'Xnk) 


n  -6x  ,  (x-x  ,  )  . 

+  l  nk  ,  ^  *4,.(x) 


k=l  1-x 


nk 


nk 


n  (29x2  -4)  (x-x  .  )2  . 

+  I  - nk  „  „  ■"*—  Jt4,  (x) 


2  2 

k=l  2(l-x.) 

nk 


nk 

3 


k=l 


2  3 

2(1-xnk) 


(-  I2> 


n  (15x  .  -35x  .  ) (x-x  ,  )  , 

V  nk  nk  nk'  „4  ,  N  ,  v  n 

+  l  - 3“ 3 -  **60  (=  U) 


(x)  , 


■ 


■ 

. 


n  -lONx  (x-x  ) 

+  I  - — 2— 2— 

k=l  3(l-x.)  n 

nk 


<■  D 


Now, 


IV  i  iv  ^nk£nk(x)^(x"Xnk)  2 

iix i  =  i  i  -njc-nk  2 — ^  &« 

k=l  1-x  , 

nk 


=  o(i) 

m 


by  remarks  (3)  and  (4) . 


II. 


n  (29x2  -4)U2(x)&2  (x) 
£  nk  n  nk 

k=l  2 (n+1) 2 


by  remarks  (2)  and  (3) . 


II. 


n  (15x  ,-35x3  )U2(x)(x-x  )  9 

IV  nk  nk  n  nk  „2  ,  N  i 

I  L  - 2 - 2 -  £nk(x)  1 

k=l  2 (n+1)  (1-x,)  nK 

nk 


=  0  (n 


-(5/2) 


by  remarks  (2)  and  (4) . 


II, 


n  lONx  ,  U2(x) (x-x  ,  ) (1-x2  )  „ 

V  nk  n  nk  nk  nl  ,  N | 

l  - 3 - 3 -  *nk(x)  I 


k=l  3  (n+1) 2  (1-x2,  ) 

nk 


since 


N 


(n+1) 


<  1,  and  also  by  remarks  (2)  and  (4).  Finally, 


2N(x-x  )2  ,  2N  U2(x) (1-x2  )£2  (x) 

nk  04  ,  >.  n  nk  nk 

- 1”  £nk(x)  =  - TT - 


3(1'Xnk) 


3 (n+1) 


so  the  lemma  is  established. 


Lemma  5.3: 


n  2N  U2  (x)  (1-x2  )£2  (x)  9  n /9s 

l  - H - nL_nk -  <1+  0(n-(l/3)) 

k=l  3 (n+1)  3 


x  e  I  , 
e’ 


uniformly  on 


as  n  -*■  00 . 


' 


■ 


■ 


Proof : 


a 


l  U^(x)  (1-x2  )SL2  (x) 
,  n  nk  nk 


l  U^(x)(l-x2,  )£2.  (x) 

/,  /,%  n  nknk 

|x-x  ,  | >n~ (1/3) 
nk '  — 


<-  Ix) 


l  U^(x)(l-x2  )£2  (x) 

_ /o\  n  nk  nk 

x-x  .  |  <n  K  ; 
nk1 


<-  L> 


Now, 


11,1  =  I  I 


Un(x)(1~xnk)3 

I  -(1/3)  (x-x  ,  ) 2 (n+1) 2 
x-x  ,  >n  nk 

nk 


^  •  n  •  (n+1)  2  •  n2^3  =  0  (n  (see  remarks 


preceding  Lemma  5.1). 

~  y  sin2( (n+l)arc  cos  x) (1-x2  ) 

Notice  that  U  (x)  (1-x  .  )  =  - - - - — 

n  nk  ,  2 


1-x 


1-x 


nk 


1-x2 


=  1  +  0(|x-xnk| ) , 


so 


L  -  l 


x-x  .  <n 
nk ' 


-(1/3) 


Un(x)(1-Xnk^nkW 


<  (l  +  0(n  (1/3)))  l  H2k(x) 


k=l 


<,  (l  +  0 (n  (1/3)))(l  +  0 (n  (1/2)))  =  1  +  0(n  (1/3)) 


2N  2 

(see  Lemma  5.1),  and  the  lemma  is  established  (since  - ^  <  —  ) . 


3 (n+1) 


2  3 


' 


■ 


, 


33. 


Lemma  5.4: 


l  >  I  +  °(n_(1/3)). 


k=l 


as  n  -►  °°, 


uniformly  for  x  e  I  . 


Proof :  This  result  is  an  immediate  consequence  of  Lemmas  5.2  and  5.3 


Lemma  5.5; 


n  2  1 

I  mnk^x)  =  °(t)  »  as  n 


k=l 


nk 


n 


uniformly  for  x  e  1^. 


Proof ; 


n  n  U  (x)  Too 

l  mnk(x)(x"xnk)  =  l  - 2  (1"Xnk)  \k(x) 

k=l  n  n  k=l  (n+1)2  nk  nk 


M  n 


(n+l)2  k=l  ~nk 


l  l*  (x)  =  0(-%)  , 


as  n  -*  00 , 


n 


by  Lemma  5.1. 


Proof  of  Theorem  5.2:  The  result  stated  in  Theorem  5.2  is  an  immediate 
consequence  of  Lemmas  5.4  and  5.5  and  Corollary  2.1. 


Theorem  5.3:  Let  x  - , • • • ,x  (-1  <x  <  •  •  •  <  x  ..  <  1)  be  the  zeros 
- nl  nn  nn  nl 

of  P  (x)  (n  =  1,2,***).  For  each  odd  n  >_  3  define  L  :  C[-l,l]  -> 
C[-l  +  n-(1/3\  1  -  n-(1/3>]  by: 


n 


L  f(x)  = 
n 


l  mnk(x)£(xnk) 
k=l 


n 


E  ”nk(x) 


k=l 


where 


' 


34. 


mnk(x) 


l-2x  .  x+x  ,  . 

_ ^ _ Sk  £4  (x) 

2  nkW’ 


1-x 


nk 


*nk(x)  = 


Pn« 


(x-x  )P' (x  ,  ) 
nk  n  nk 


-(1/3)  -(1/3), 

,  1  -  n  J . 


Let  11*1^  a,  denote  the  sup  norm  on  [-1  +  n 
Then  there  exists  a  constant  A  >  0  and  an  integer  n^  so  that  for 
every  odd  n  _>  n^  and  every  f  e  C[-l,l], 


L  f-f  |L  <  A  «o_  f— )  . 

n  "n,°°  —  f  hi' 


Remarks :  Notice  that  the  error  estimate  applies  on  an  increasing 

sequence  of  intervals  whose  union  is  (-1,1).  However,  that  this 
cannot  be  improved  so  that  the  error  estimate  applies  to  [-1,1]  is 
demonstrated  by  the  following: 


It  will  now  be  shown  that  there  exists  a  function  f  e  C[-l,l] 


s.t.  L  f(l)  +  f(l). 
n 

If  k  £  y  ,  then  (Szego  [10],  p.122) 


(k-^)  it 
n 


<  arc  cos  x 


kr 

nk  —  n+1 


9 


and  hence  there  exist  constants  >  0  s.t.  for  every  n,k  with 


k  if  , 


c,£)2n-xnk<Cl£)2  . 


Also  (Szegc5  [10,  p.238),  if  k£y  , 


%  k 


-(3/2)  2 
n 


9 


■ '  '■ 


35. 


and 


1  <  1  +  x  .  <  2  , 
—  nk  — 


so  if  k  <  n/2,  «*(!>  -  -■ 2,nV  Xnk  *  -T  • 

1  -  Xnk 


Thus  there  exist  constants  d  >  >  0  s.t.  for  every  n,k  with 

k-f  ’ 


dl  d2 

2  >  mnk(1)  >  2  * 

n  n 

d2  n 

Let  0  <  a  <  1  s.t.  card  {k|xnk  _>  a}  <_  — —  ,  for  every  n  =  1,2, 


Define 


1  ,  x  =  1 


f(x)  =  /  0  ,  x  e  [-l,a] 


x-a  ,  . 

ITJ-  ,  X  6  (a,l) 


Then  f  e  C [-1 , 1]  and 


and  so 


l  “T  f  (Xnk) 

11  I 

Lnf (1)  -  ~^~n  d~ - ±¥’ 

l  ~i 

k=l  n 


L  f(l)  +  f(l) 
n 


Proof  of  Theorem  5.3:  Some  of  the  simpler  details  will  be  omitted  from 
the  proof. 

First,  let  E,  <•••<£  be  the  zeros  of  P'(x). 

n,n-l  nl  n 

It  is  easily  proven  that,  for  each  n  >  3  and  2  _<  k  _<  n-1. 


(I) 


mm 


XeI-?nk,?n,k-l] 


£  .  (x)  =  min{ £  .  (5  .  ),  Z  ,  (?  ,  ,)}  >  0, 

nk  nk  nk  nk  n,k-l 


. 
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Certain  inequalities  will  be  needed,  and  are  listed  below: 

(1)  (Balazs  and  Turan  [2]).  If  n  is  odd,  n  3,  then 

lPn^nk^l  -  (8?rJnk)  ( 1/2\  where  =  min  (k,n+l-k) . 

(2)  (Szego  [10],  p.238) . 


|P'(x  k)|  ^  j  n2 

n  nk  1  Jn+l,k 


(3)  (Szego  [10],  p.122). 


(k-^ir  „  kTT  rn-, 

-  <  arc  cos  x  .  <  — —  ,  k  <  [— ] , 

n  —  nk  —  n+1  —  2 


and  a  symmetric  result  holds  for  the  values  k  >  y  . 


Easy  calculations  yield 


Xnk  -  Vk+l1  -  (4jnk+1)  h  (k  = 

n 


and 


1  -  Xnk  i 


cA  \  2 
r3n+l,k 


n 


(k  =  1, * • • ,n) . 


(4)  (Szego  [10,  p.165). 


Pn(x)I  i 


Jf  n-(1/2>(l-x2)-(1/4) 


Using  these  inequalities,  it  can  be  shown  that  there  exists 
a  constant  C  >  0  s.t.,  for  every  odd  n  3, 


(id 


mlnUnk(5nj) 


k  =  2,***,n-l;  j  =  k-l,k}  _>  C. 


For  example,  let  2  <_  k  <_  n-1.  Then  there  exists  a  constant  >  0  s.t 


■■ 


37. 


‘nk^nk3  ‘ 


P  (5  .  ) 
n  nk 


£  , -x  P'  (x  ) 
nk  nk  1  n  nk 


>  C 


°  ,,a  ,1X  2  2 

(4jnk+1)  77  n 


2/2  7T5/2 


,3/2 
J  n+1 , k 


o 


'(A3nk+1)/V 


40tt 


5/2 


>  0. 


The  treatment  of  i  ,  (£  .  .)  is  similar. 

nk  n,k-l 

From  (I)  and  (II)  it  follows  by  easy  calculations  that  there 
exists  >  0  s.t.  for  odd  n  >_  3, 


(HI) 


n 


min  /  m  ,  (x)  >  C-  . 

rr  r  1  1  1  nk  ! 

xc  U  ,  ,£  , ]  k=l 
n,n-l  nl 


Finally,  we  will  show  that 


n 


2  1 

(IV)  sup  _(1/3)  l  mnk(x)  (x_xnk)  =  ’  odd 

x  <l-n  v  '  k=l  n 


n 


Let  n  be  odd,  n  >  3,  and  x  e  [— 1+n  ^  ^  ^,1-n  ^  ^  The  sum  will 


n 


-(1/3) 


be  considered  in  two  parts  -  those  k  for  which  x-x  .  <  - » -  , 

n-(1/3)  .  .  ..  1  -(1/3) 


and  those  k  for  which  x-x  .  > 

1  nk 1 


.  Let  a  =  —  n 
n  2 


If 


hnk(x)  = 


1  -  2x  ,  x  +  x  . 

nk  nk 


1-x 


^nk^x3  ’ 


nk 


then 


n 


l  hnk(K)  E 


k=l 


I  mnk(x)(x-xmi,)2  =  .  I 


h  (x)P  (x) 
nk  n 


x-x  ,  <a 
nk n 


nk' 


c-x  ,  |  <a  (p*  (x  )) 
nk1—  n  v  n  nk  ' 


0(4) 


V  1  ,  hnkW 

n  x-x  .  <a 
1  nk1—  n 


^n+1 ,k 


n 


a-x2)-(1/2) 


■ 


■ 


* 
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=  0("t).  I.  hnk(x) 


'3n+l,kl 


n  x-x  ,  <a 
nk1—  n 


n 


(since  (1-x2)  (1/2)  =  (1-x2  )  (1/2) 

nk 


r,  2 

1-x 


1-x 


(1/2) 


=  o(a^k)-(1/2>).) 


nk 


=  hnk(x) 


n  k=l 


■  • 
n 


Also,  quite  easily  one  gets 


I 


x-x  .  >a 
nk 1  n 


mnk(x)(K-xnk) 


(l-2x  x+x2  JP^Cx) 
nk  nk  n 


l 

X~Xnk  I  >an  (1-<k)  (5Mcnk)  ^  (Pn  (xnk) ) 


=  0(4)  - 


n 


establishing  IV. 


Notice  that  certainly 


c  [E  ..5.], 
n,n-l  nl 


so  Theorem  5.3  now  follows  immediately  from  III,  IV,  and  Corollary  2.1. 


■ 

. 


CHAPTER  6 


CONCLUSION 


As  predicted  in  Chapter  1,  the  error  estimates  achieved  in 

Chapter  5  are  an  improvement  over  those  obtained  for  H-F  interpolation. 

However,  it  is  not  known  whether  the  results  of  Chapter  5  can  be 

improved.  It  is  a  simple  matter  to  show  that  if  {Lr}  is  a  sequence  of 

proper  interpolation  operators  defined  on  C [ —1 , 1 ]  and  {e  }  is  a 

n 

sequence  s.t.  e  >0,  e  10,  with 

n  n  ’ 


lim  sup 
n-*50  feC[-l,l] 


|| L  f-f 
"  n 


a)  (e  ) 

f  n 


=  0(1) 


9 


then  e  >  K  •  —  for  some  constant  K  and  all  n.  (U  •  denotes  the 
n  —  n  Nl  11  oo 

sup  norm  on  any  nontrivial  sub interval  I  of  [-1,1].).  However,  it  is 

not  known  whether  e  ^  —  is  the  best  possible  in  order  that 

n  n 


sup  lim 
feC[-l,l]  n-^°° 


L  f-f  || 
n  11 

03  (e  ) 

t  n 


=  0(1) 


even  if  {L^}  is  restricted  to  some  smaller  class  of  interpolation 
operators  which  contains  those  of  Chapter  5. 

It  seems  likely  that  the  proof  of  Lemma  4.1  could  be 
strengthened  so  that  the  /£n  £n  n  factor  could  be  eliminated.  It 
also  seems  likely  that  all  the  results  of  Chapter  5  could  be  combined 
into  a  general  theorem  using  ultraspherical  polynomials  as  the  building 
blocks  for  some  rational  interpolation  operator. 
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